Accepted by ApJ: February 17, 2010 



Local Benchmarks for the Evolution of Major-Merger Galaxies — 
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ABSTRACT 

We present Spitzer observations for a sample of close major- merger galaxy 
pairs (KPAIR sample) selected from 2MASS/SDSS-DR3 cross-matches. The 
goals are to study the star formation activity in these galaxies and to set a 
local bench mark for the cosmic evolution of close major mergers. The Spitzer 
KPAIR sample (27 pairs, 54 galaxies) includes all spectroscopically confirmed 
spiral-spiral (S+S) and spiral-elliptical (S+E) pairs in a parent sample that is 
complete for primaries brighter than K=12.5 mag, projected separations of 5 < 
s < 20h~ 1 kpc, and mass ratios < 2.5. The Spitzer data, consisting of images 
in 7 bands (3.6, 4.5, 5.8, 8, 24, 70, 160/xm), show very diversified IR emission 
properties. Compared to single spiral galaxies in a control sample, only spiral 
galaxies in S+S pairs show significantly enhanced specific star formation rate 
(sSFR=SFR/M), whereas spiral galaxies in S+E pairs do not. Furthermore, 
the SFR enhancement of spiral galaxies in S+S pairs is highly mass-dependent. 
Only those with M £ 1O 1O ' 5 M show significant enhancement. Relatively low 
mass (M ~ 1O 1O M ) spirals in S+S pairs have about the same SFR/M compared 



1 Infrared Processing and Analysis Center, California Institute of Technology 100-22, Pasadena, CA 91125 

2 Gcorgia College & State University, CBX 82, Milledgeville, GA 31061 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

4 Purple Mountain Observatory, Chinese Academy of Sciences, 2 West Beijing Road, Nanjing 210008, 
China 

institute of Astronomy, National Central University, Chung-Li 32054, Taiwan 

institute of Astrophysics and Leung Center for Cosmology and Particle Astrophysics, National Taiwan 
University, Taipei, Taiwan 10617 

7 Spitzer Science Center, California Institute of Technology, Mail Stop 220-6, Pasadena, CA 91125 



-2 - 



to their counterparts in the control sample, while those with lO n M have on 
average a ~ 3 times higher SFR/M than single spirals. There is evidence for 
a correlation between the global star formation activities (but not the nuclear 
activities) of the component galaxies in massive S+S major-merger pairs (the 
"Holmberg effect"). There is no significant difference in the SFR/M between 
the primaries and the secondaries, nor between spirals of SEP < 1 and those of 
SEP > 1, SEP being the normalized separation parameter. The contribution of 
KPAIR galaxies to the cosmic SFR density in the local universe is only 1.7%, 
and amounts to PkPAiR = 2 -54 x 10" 4 (M yr" 1 Mpc~ 3 ). 

Subject headings: galaxies: interactions — galaxies: evolution — galaxies: star- 
burst — galaxies: general 

1. Introduction 

Two questions about the evolution of galaxy mergers are being intensely debated in 
the current literature: (1) Does merger rate have a strong or weak cosmic evolution? (2) 
Are mergers responsible for the strong evolution of the cosmic star formation rate (SFR) 
since z ~ 1 (e.g., Lilly et al. 1996; Madau et al. 1998; Hopkins 2004)? Answers to these 
questions have important implications to the understanding of basic processes in galaxy 
formation/evolution, such as mass growth and star formation. 

Earlier studies of merger rate, using samples of close galaxy pairs and morphologically 
disturbed systems in different redshift ranges, yielded a very broad range of the evolutionary 
index m, m — - 6, assuming the evolution has a power-law form (l + z) m (Zepf & Koo 1989; 
Burkey et al. 1994; Carlberg et al. 1994; Yee & Ellington 1995; Woods et al. 1995; Patton 
et al. 1997; Wu & Keel 1998; Brinchmann et al. 1998; Le Fevre et al. 2000; Carlberg et al. 
2000; Patton et al. 2002; Conselice et al. 2003; Lin et al. 2004; Bundy et al. 2004; Lavery et 
al. 2004). More recent results can be divided into two camps, the "strong evolution" camp 
with m ~ 3 (Conselice 2006; Kampczyk et al. 2007; Kartaltepe et al. 2007) and the "weak 
evolution" camp with m ~ 0.5 (Lin et al. 2008; Lotz et al. 2008). The strong evolution 
scenario is consistent with the cosmic time dependence of the merging rate of dark matter 
halos (Lacey & Cole 1993; Khochfar & Burkert 2001). However, more recent simulations 
including sub-halo structures support the weak evolution scenario (Berrier et al. 2006; Guo 
& White 2008). 

Since the discovery of a strong evolution of the cosmic SFR (Lilly et al. 1996; Madau 
et al. 1998), many authors have argued that the primary cause is a rapid decline of merger- 



-3 - 



induced star formation (Driver et al. 1995; Glazebrook et al. 1995; Abraham et al. 1996; 
Brinchmann et al. 1998, Le Fevre et al. 2000; Conselice et al. 2003). In particular, infrared 
(IR) surveys by both ISO (Elbaz et al. 2002) and Spitzer (Le Floch et al. 2006) found 
that beyond z ~ 0.7-1.0, the cosmic star formation is contributed mostly by luminous IR 
galaxies (LIRGs) with Lm > 1O 11 L . In the local universe most LIRGs are in merger 
systems, and they contribute only a few percent of the integrated IR luminosity density 
(Sanders & Mirabel 1996). This seems to provide another argument for a strong evolution 
of merger-induced starbursts as the primary driver of the strong cosmic SFR evolution. 
However, are LIRGs at z ~ 1 indeed mostly in merger systems, as observed for their local 
counterparts? To this question there are both positive answers (Zheng et al. 2004; Hammer 
et al. 2005; Bridge et al. 2007) and negative answers (Flores et al. 1999; Bell et al. 2005; 
Melbourne et al. 2005; Lotz et al. 2008). Authors in the latter group found that most 
LIRGs at z ~ 1 are normal late-type galaxies. According to them, it is the secular evolution 
of normal late-type galaxies that is mostly responsible for the cosmic SFR evolution, not the 
evolution of mergers. 

What are the reasons for these controversies? The foremost among them are sample 
selection effects. There are two classical methods of selecting interacting/merging galaxies. 
One is to find galaxies with peculiar morphology (e.g. with tidal tails, double nuclei, or 
distorted discs), and the other is to identify paired galaxies. There are systematic differences 
between interacting galaxies selected using these two different methods. In the so-called 
merger sequence (Toomre 1977), galaxies in close pairs are usually mergers in the early 
stages when the two galaxies are still separable. In contrast, peculiar galaxies are mostly 
found in the later stages when the first collision between the two galaxies is happening or 
passed. However, this distinction is not clear-cut. Mergers such as the Antennae (= Arp 244) 
can be identified using both methods. 

There are several known biases in the morphological selection method. The most seri- 
ous one is due to the mis-identifications of isolated irregular galaxies or starburst galaxies 
as mergers (Lotz et al. 2008). This effect is particularly severe for samples selected in 
the rest frame UV bands (and to some extend those selected in the rest frame B-band), 
where the emission of young stars and dust extinction can significantly affect the surface 
brightness distribution. Another bias is caused by missing low surface brightness merger 
features, such as faint tidal tails and bridges, in observations that lack sufficient sensitivity. 
This incompleteness becomes increasingly severe for high-z surveys because of the cosmic 
dimming. 

The most common bias affecting current pair selected samples is an incompleteness 
known as "missing the secondary" (Xu et al. 2004). For flux limited (= apparent magnitude 
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limited) samples or luminosity limited (= absolute magnitude limited = "volume limited") 
samples, a paired galaxy brighter than the limit can be missed if its companion is fainter 
than the limit. The amplitude of this incompleteness can vary with the redshift, and cause 
significant bias in the results on the evolution of merger rate. For example, in many recent 
studies of merger rate evolution, pair fractions of galaxies of M < Mu m are compared, where 
Mi im is an absolute magnitude limit in the rest frame B or V band. For these samples, the 
amplitude of the incompleteness due to the "missing the secondary" bias is on the order of 
Q ~ 0.5<f){M lim )5M/ J^J 4>(M) dM, where is the luminosity function (LF) of the paired 
galaxies (e.g. Xu et al. 2004; Domingue et al. 2009), and 5M is the typical magnitude 
difference between the two galaxies in a pair, which is ~ 1 mag for major mergers. When 
Mn m is fixed, the ratio Q decreases with z if the LF has a positive "luminosity evolution" 
(i.e. the "knee" of the LF becomes brighter with increasing z), as being observed for the 
LF of field galaxies (Wolf et al. 2003; Marshesini et al. 2007). This can introduce an 
artificial "evolution" of the merger rate in studies comparing pair fractions in samples at 
different redshifts that are limited by the same Mn m (e.g. Kartaltepe et al. 2007), and being 
responsible for the high evolutionary index (m~ 3) found in those studies. On the other 
hand, all studies in the 'weak evolution' camp invoke a correction for the "passive luminosity 
evolution" (PLE), allowing the absolute value of Mu m to increase with z accordingly. This 
indeed reduces the effect of the incompleteness on the evolutionary index of the merger rate. 
However, as pointed out by Kartaltepe et al. (2007), there is no strong empirical justification 
for the PLE model. If the true evolution of the luminosities of interacting galaxies is different 
from PLE, then the bias is still present. It is better to get rid of the incompleteness from 
the merger rate studies in the first place. 

Other biases for pair selections include: (1) Contamination due to unphysical, projected 
pairs ("interlopers"). This affects mostly samples with incomplete redshifts or with only 
photometric redshifts. (2) Incompleteness due to missing of pairs in which the two galaxies 
are too close to be separated visually because of insufficient angular resolution or obscuration 
by dust. 

Being aware of the selection effects that lead to the conflicting results, we set out to 
design a set of merger selection criteria that minimize the biases mentioned above. Firstly, we 
opted for the pair selection method instead of the morphological selection method because, 
based purely on galaxy separation, pair selections are more objective. However, this also 
confines our study to systems prior to the final stages of the merging process. Secondly, 
we chose the rest frame K as the waveband in which our samples are to be selected. This 
is the band least affected by star formation and dust extinction and most closely related 
to mass (Bell & De Jong 2001). Thirdly, we confine ourselves to close major-merger pairs 
with mass ratio < 2.5 and with projected separation in the range 5h _1 < s < 20h _1 kpc 
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(h = H /(100 km sec _1 Mpc -1 )). Many studies have found that only major mergers with 
separations comparable to the size of galaxies (i.e. £ 20h~ 1 kpc) show significant SFR 
enhancements (Xu & Sulentic 1991; Barton et al. 2000; Lambas et al. 2003; Nikolic et al. 
2004; Alonso et al. 2004; Woods et al. 2006; Barton et al. 2007; Ellison et al. 2008). Other 
detailed selection criteria are presented in Section 2. It should be emphasized that, when 
studying the evolution of interacting galaxies, samples at different redshifts must be selected 
using identical criteria. 

In two earlier papers (Xu et al. 2004; Domingue et al. 2009), we started from samples 
in the local universe, with the goal of setting local benchmarks for merger rate evolution. 
Xu et al. (2004), using a sample of 19 close major-merger pairs selected from the matched 
2MASS/2dFGRS catalog (Cole et al. 2001), derived the K-band luminosity function (KLF) 
and the differential pair fraction function (DPFF) of local binary galaxies. This was followed 
up by a more extended analysis (Domingue et al. 2009), exploiting a large sample of 173 
close major- merger pairs selected from 2MASS/SDSS-DR6 cross-matches. Assuming the 
mass dependent merger time scale of Kitzbichler & White (2008), Domingue et al. (2009) 
found that the differential merger rate increases with mass, and the merger rate v.s. mass 
relation is in good agreement with what being found in the N-body simulations of Mailer et 
al. (2006). 

In this paper we report Spitzer imaging observations (7 bands at 3.6, 4.5, 5.8, 8, 24, 
70, and 160yum) of 27 galaxy pairs in the local universe, selected from cross matches of 
2MASS and SDSS-DR3. The scientific goals are (1) studying the star formation activity 
in these galaxies and (2) setting a local bench mark for the evolution of the SFR in close 
major mergers. Assuming that only late type galaxies contribute significantly to the total 
SFR, our analysis is concentrated on spiral (S) galaxies in spiral-spiral (S+S) and spiral- 
elliptical (S+E) pairs. The main focus of this paper is on the enhancement (or the lack 
of it) of the SFR of paired galaxies. Most previous studies on this subject are based on 
optical/UV observations susceptible to dust extinction (Barton et al. 2000; Lambas et al. 
2003; Nikolic et al. 2004; Alonso et al. 2004; Woods et al. 2006; Woods & Geller 2007; 
Barton et al. 2007; Ellison et al. 2008). Early FIR studies based on IRAS observations 
(Kennicutt et al. 1987; Telesco et al. 1988; Xu & Sulentic 1991) can not resolve the pairs 
because of the coarse angular resolution of IRAS. The more recent ISO study of Xu et al. 
(2001) and Spitzer study of Smith et al. (2007), which resolved pairs into discrete regions 
(e.g. nuclei of component galaxies, bridges and tails, overlapping regions between the two 
disks, etc.), are confined to pairs with strong interacting features, and therefore are biased to 
interacting galaxies with strong star formation enhancement. This work shall be neutral to 
these biases. Our studies on IR SED's (Domingue et al. in preparation) and on the optical 
properties including new Ha and H/3 imaging observations (Cheng et al. in preparation) will 
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be published separately. Throughout this paper, we adopt the A-cosmology with Q m = 0.3 
and tt A = 0.7, and H = 75 (km sec^Mpc -1 ). 

2. The Sample 

The sample of galaxy pairs was selected from cross-matches between the Two Micron 
All Sky Survey (2MASS) Extended Source Catalog (XSC; Jarrett et al. 2000) and the galaxy 
catalog of Sloan Digital Sky Survey (SDSS) Data Release 3 (DR3; Abazajian et al. 2005). 
The selection procedure is similar to that in Xu et al. (2004) and Domingue et al. (2009). 
First, the following two criteria were set for galaxies to be considered in the pair selection: 

(1) Galaxies should be brighter than i^ s =13.5 mag, the completeness limit of the XSC 
(Jarrett et al. 2000). The default K20 magnitude is used for the K s band (2.16/zm) 
fluxes (Jarrett et al. 2000). The K20 magnitudes were taken from Domingue et al. 
(2009), for which the photometric error due to the blending of close neighbors was 
corrected. As in Xu et al. (2004), a uniform -0.2 mag correction was applied to the 
K 20 of galaxies when extrapolation to total K band magnitude was necessary. 

(2) Galaxies should have the redshift completeness index c 2 >0.5, where c 2 is the ratio of 
the number of galaxies with measured redshifts within 1° radius from the center of 
the galaxy in question and the number of all galaxies within the same radius (Xu et 
al. 2004). This confines the pair selection to regions where the SDSS-DR3 has good 
spectroscopic coverage. 

The resulting sample has 50312 galaxies, of which 42847 have measured redshifts (85% 
redshift completeness). A comparison with the number counts of 2MASS galaxies (Kochanek, 
et al. 2001) yielded an equivelent sky coverage of Q = 3000 deg 2 (Q = N/CN, where number 
counts CN is in units of deg -2 ). 

Pair selection was then carried out with the following criteria: 

(1) The K s magnitude of the primary is brighter than 12.5 mag. A primary is defined as 

the brighter component of a pair. 

(2) The K s difference between the two components is less than 1 magnitude: <5K S < 1 mag. 

(3) At least one of the components has spectroscopic redshift. 

(4) The projected separation is in the range of 5 < s < 20h _1 kpc. 



-7- 



(5) The velocity difference is less than 500 km sec 1 : 5v < 500 km sec 1 . 

Criteria (1) and (2) ensure that, by construction, all galaxies in this pair sample are brighter 
than K s =13.5, the completeness limit of 2MASS survey (Jarrett et al. 2000). Therefore 
the sample does not suffer from the "missing the secondary" bias that plagued many earlier 
merger studies (see Section 1). This also restricts our sample to major merger pairs whose 
mass ratios are < 2.5. Criterion (4) makes the comparisons between local pairs in this sample 
and pairs in high-z samples robust. For high-z galaxies, it is difficult to distinguish a galaxy 
pair from a single galaxy if the separation is < 5h~ l kpc, and the probability for chance pairs 
is significantly higher if projected separation is much larger than 20h~ 1 kpc (Kartaltepe et 
al. 2007). 

A total of 57 pairs were selected according to the criteria. Component galaxies are 
classified as either "S" or "E" using the following scheme: First, classical "eye-ball" clas- 
sifications were carried out. Independently, two of us (CKX and DD) inspected the SDSS 
optical image of every galaxy in the sample and assigned a type to it according to its mor- 
phology. In addtion, we also run an automatic classification script. It classifies a galaxy 
as "E" if it — r > 2.22 and R 5 o/R 9 q < 0.35, where i? 50 and R 90 are the radii containing 
50% and 90% of the Petrosian flux (Park & Choi 2005). Otherwise, the galaxy is classified 
as "S" type. The final classification is the median of the two eyeball results and the result 
of the automatic method. We excluded 20 pairs in which both components are "E" types 
("E+E" pairs). According to their optical colors and EW(Ha), little star formation is evi- 
dent, consistent with IRAS results on E+E pairs (Xu & Sulentic 1991). In order to minimize 
the contamination due to interlopers, single-redshift pairs were also excluded. We would 
stress that, judging from SDSS images, these single-redshift pairs do not show any special 
characteristics compared to the total sample, therefore dropping them will not introduce any 
significant bias to our study. 

The final sample (KPAIR hereafter) for the Spitzer observations contains 27 S+S and 
S+E pairs. Among the 54 component galaxies, 42 are classified as S and 12 as E. They 
are listed in Table 1. The mass in the table (and hearafter) is the so called "stellar mass", 
excluding the mass of the dark matter and gas. It is estimated from the K-band luminosity 
by assuming a mass-to-luminosity ratio of 1.32 (solar units), taken from Cole et al. (2001) 
for the Salpeter IMF[j] The heliocentric velocities are from SDSS-DR3, and the 60/xm flux 
densities are from the IRAS Faint Source Catelog (Moshir et al. 1992). 



1 The mass-to-luminosity ratio estimated using the Kennicutt IMF is a factor of 0.55 lower (Cole et al. 
2001), and that estimated using the Kroupa IMF is a factor of 0.59 lower (Kauffmann et al. 2003). 
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3. Spitzer Observations 



Infrared imaging observations of all pairs but J1315+6207 in Table 1 were carried out in 
2005 and 2006 under the Spitzer Cycle 2 GO Program "Local Benchmarks for the Evolution 
of Interacting Galaxies" (PID #20187). These include images in the four bands (3.6, 4.5, 
5.8 and 8.0 /an) of the Infrared Array Camera (IRAC, Fazio et al. 2004) and in three bands 
(24, 70 and 160 fim) of the Multiband Imaging Photometer for Spitzer (MIPS, Rieke et 
al. 2004). KPAIR J1315+6207 (UGC08335a/b) is included in the Spitzer GO-1 survey of 
LIRGs (Mazzarella et al. 2009). Its IRAC data were taken from Mazzarella et al. (2009), 
while its MIPS observations were unsuccessful. 
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Table I. KPAIR Galaxy Sample 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Pair ID 


Galaxy ID 


RA 


Dec 


V z 


K s 


log(M) 


Type 






(KPAIR) 


(2MASX) 


(J2000) 


(J2000) 


(km s _1 ) 


(mag) 


(Mq) 




(Jy) 


(log(L©)) 


J0020+0049 


J00202580+0049350 


00h20m25.8s 


+00d49m35s 


5078 


10.99 


10.84 


S 


0.62 


9.68 




J00202748+0050009 


00h20m27.5s 


+00d50m01s 


5480 


10.50 


11.04 


E 






J0109+0020 


J01093371+0020322 


01h09m33.7s 


+00d20m32s 


13499 


12.39 


11.08 


E 








J01093517+0020132 


01h09m35.1s 


+00d20ml3s 


13319 


12.47 


11.05 


S 






J0118-0013 


J01183417-0013416 


01hl8m34.1s 


-00dl3m42s 


14219 


12.05 


11.27 


S 2 


3.48 


11.30 




J01 183556-0013594 


01hl8m35.6s 


-00dl3m59s 


14273 


12.88 


10.93 


S 






J0211-0039 


J02110638-0039191 


02hllm06.3s 


-00d39m21s 


5920 


11.42 


10.77 


S 2 








J02110832-0039171 


02hllm08.3s 


-00d39ml7s 


6016 


10.90 


10.98 


S 






J0906+5144 


J09060283+5144411 


09h06m02.8s 


+51d44m41s 


8849 


11.68 


11.01 


E 








J09060498+5144071 


09h06m05.0s 


+51d44m07s 


8852 


11.95 


10.90 


S 2 






J0937+0245 


J09374413+0245394 


09h37m44.1s 


+02d45m39s 


6988 


10.01 


11.46 


S 


2.00 


10.43 




J09374506+0244504 


09h37m45.0s 


+02d44m50s 


6790 


10.45 


11.29 


E 






J0949+0037 


J09494143+0037163 


09h49m41.4s 


+00d37ml6s 


1861 


11.59 


9.71 


S 


2.27 


9.36 




J09495263+0037043 


09h49m52.6s 


+00d37m05s 


1918 


10.98 


9.95 


s 






J1020+4831 


J10205188+4831096 


10h20m51.9s 


+48d31ml0s 


15886 


13.26 


10.88 


s 








J10205369+4831246 


10h20m53.6s 


+48d31m24s 


15930 


12.27 


11.27 


E R 






J1027+0114 


J10272950+0114490 


10h27m29.5s 


+01dl4m48s 


6727 


11.79 


10.73 


s 








J10272970+0115170 


10h27m29.7s 


+01dl5ml6s 


6661 


10.90 


11.08 


E 






J1043+0645 


J10435053+0645466 


10h43m50.5s 


+06d45m47s 


8262 


11.96 


10.83 


S 








J10435268+0645256 


10h43m52.7s 


+06d45m25s 


8088 


12.20 


10.73 


s 






J1051+5101 


J10514368+5101195 


10h51m43.6s 


+51d01m20s 


7503 


10.27 


11.41 


E 


0.78 


10.07 




J10514450+5101303 


10h51m44.5s 


+51d01m30s 


7138 


10.97 


11.13 


S 






J1202+5342 


J12020424+5342317 


12h02m04.3s 


+53d42m32s 


19366 


12.97 


11.16 


s 








J12020537+5342487 


12h02m05.3s 


+53d42m48s 


19156 


12.43 


11.37 


E 






J1308+0422 


J13082737+0422125 


13h08m27.4s 


+04d22ml3s 


7186 


13.39 


10.15 


S 








J13082964+0422045 


13h08m29.6s 


+04d22m05s 


7251 


12.44 


10.53 


s 






J1332-0301 


J13325525-0301347 


13h32m55.3s 


-03d01m35s 


14297 


12.95 


10.90 


s 


0.57 


10.51 




J13325655-0301395 


13h32m56.6s 


-03d01m40s 


14000 


12.19 


11.21 


s 






J1346-0325 


J13462001-0325407 


13h46m20.0s 


-03d25m41s 


6949 


11.20 


11.01 


s 








J13462215-0325057 


13h46m22.2s 


-03d25m06s 


7171 


11.66 


10.82 


E 2 






J1400+4251 


J14005782+4251207 


14h00m57.8s 


+42d51m20s 


9689 


11.89 


11.01 


S 


2.32 


10.80 




J14005882+4250427 


14h00m58.8s 


+42d50m42s 


9923 


12.18 


10.90 


S 






J1425+0313 


J14250552+0313590 


14h25m05.5s 


+03dl3m59s 


10693 


11.98 


11.06 


E 1 








J14250739+0313560 


14h25m07.4s 


+03dl3m55s 


10807 


12.97 


10.66 


S 






J1433+4004 


J14334683+4004512 


14h33m46.8s 


+40d04m52s 


7674 


10.78 


11.25 


S 


1.80 


10.48 




J14334840+4005392 


14h33m48.4s 


+40d05m39s 


7788 


11.17 


11.10 


s 






J1453+0317 


J14530282+0317451 


14h53m02.8s 


+03dl7m46s 


1576 


10.67 


9.92 


s 


0.29 


8.30 




J14530523+0319541 


14h53m05.2s 


+03dl9m54s 


1573 


10.05 


10.17 


s 






J1506+0346 


J15064391+0346364 


15h06m43.9s 


+03d46m36s 


10498 


11.48 


11.22 


s 


0.27 


9.91 




J15064579+0346214 


15h06m45.8s 


+03d46m22s 


10183 


11.61 


11.17 


s 






J1510+5810 


J15101587+5810425 


15hl0ml5.8s 


+58dl0m43s 


9161 


11.77 


11.02 


s 


0.29 


9.86 




J15101776+5810375 


15hl0ml7.8s 


+58dl0m37s 


9563 


12.35 


10.79 


s 






J1528+4255 


J15281276+4255474 


15h28ml2.8s 


+42d55m48s 


5588 


10.02 


11.26 


s 


1.10 


9.96 
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Table 1 — Continued 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Pair ID 


Galaxy ID 


RA 


Dec 




K s 


log(M) 


Type 






(KPAIR) 


(2MASX) 


(J2000) 


(J2000) 


(km s _1 ) 


(mag) 


(M Q ) 




(Jy) 


(log(L )) 




J15281667+4256384 


15h28ml6.6s 


+42d56m39s 


5345 


10.59 


11.03 


S 






J1556+4757 


J15562191+4757172 


15h56m21.9s 


+47d57ml7s 


5741 


12.10 


10.49 


S 


0.53 


9.71 




J15562738+4757302 


15h56m27.4s 


+47d57m30s 


5980 


12.16 


10.46 


E 






J1602+4111 


J16024254+4111499 


16h02m42.6s 


+41dllm50s 


10019 


11.69 


11.11 


S 


1.43 


10.60 




J16024475+4111589 


16h02m44.7s 


+41dllm59s 


9950 


12.50 


10.78 


S 






J1704+3448 


J 17045089+3448530 


17h04m50.9s 


+34d48m53s 


17163 


13.07 


11.01 


s 


1.41 


11.06 




J17045097+3449020 


17h04m50.9s 


+34d49m02s 


16893 


12.40 


11.28 


s 






J2047+0019 


J20471908+0019150 


20h47ml9.1s 


+00dl9ml5s 


4209 


9.08 


11.37 


s 


0.34 


9.18 




J20472428+0018030 


20h47m24.3s 


+00dl8m03s 


3795 


9.74 


11.10 


E 






J1315+6207 


J13153076+6207447 


13hl5m30.8s 


+62d07m45s 


9170 


11.99 


10.91 


S 


11.39 


11.47 




J13153506+6207287 


13hl5m35.1s 


+62d07m29s 


9176 


11.54 


11.09 


S 







Note. — Descriptions of Columns: 

(1) Pair ID. The designations are "KPAIR J0020+0049", etc. 

(2) Galaxy ID, taken from 2MASS. 

(3) RA (J2000). 

(4) Dec (J2000). 

(3) Recession velocity taken from SDSS. 

(6) K s (K 2 o) magnitude taken from 2MASS. 

(7) Logarithm of the stellar mass of the galaxy. 

(8) Morphological type. Galaxies containing known AGNs (via NED) are flagged with 1, 2 or R to indicate type 1, type 2 or 

radio AGNs. 

(9) The IRAS 60/im flux density of the total pair. 

(10) The 60/jto luminosity (uh^) of the total pair. 
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3.1. IRAC 

3.1.1. Observations 

All but 3 pairs in our sample are smaller than 2'. 5 at optical wavelengths. Each pair 
(or each pair component for the 3 large pairs) was observed with the IRAC instrument in 
Full Array mode. High dynamic range settings and a 12 positions Reuleaux dither pattern 
were utilized. The dynamic range setting was designed to incorporate short exposures of 
0.4 sec along with longer exposures of 12 sec, in order to acquire photometry on potentially 
saturated nuclei as well as reach the requested signal-to-noise ratio for the outer galaxy 
regions. The dither pattern produces the redundancy needed to eliminate cosmic rays. The 
average surface brightness sensitivity levels (4<r) are 0.032, 0.048, 0.160, and 0.316 MJy sr _1 
in the 4 IRAC bands, respectively. 

3.1.2. Data Reduction 

The Basic Calibrated Data (BCD) images were produced in the Spitzer Science Center 
(SSC) pipeline and these served as the start of our data reduction process. Custom IDL 
tools created by one of us (J.S. Huang of the IRAC instrument team) were used to apply the 
additional steps of pointing refinement, distortion correction, and mosaicking with sigma- 
clipping for the rejection of cosmic rays. The long and short exposures were processed 
independently for the purpose of creating mosaics in all four bands. These mosaics were 
resampled from the native 1'.'22 pixel scale to a 0Y66 pixel scale. 

Aperture photometry was carried out on the images in order to obtain flux densities in 
the four IRAC bands. Two different methods were applied to different sources: 

(1) Standard aperture photometry for individual galaxies. This was performed on well 
separated pairs, utilizing both IPAC Skyview and IRAF APPHOT. Several blank sky 
regions were chosen from within the IRAC field containing the galaxy pairs. These 
regions were used to calculate the mean background value for subtraction. The stan- 
dard deviation of their means provides estimate of the uncertainties of background 
subtraction (Smith et al. 2007). Apertures of different shapes (circular, eliptical and 
polygon) are chosen according to the situation. 

(2) Aperture photometry first for the total pair and then for one of the components. This 

was performed on pairs whose two component galaxies are too close to be separated 
cleanly. We always chose the galaxy which has a more regular morphology (point- 
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source-like or a regular elliptical) for the aperture photometry. This flux is then sub- 
tracted from the flux of the total pair, which yields the flux of the other component. 

In addition to total fluxes of whole galaxies, circular aperture fluxes on two physical 
scales corresponding to D = 4 and 10 kpc, respectively, were also obtianed. These fluxes 
sample the nucleus and disk emissions. 

In the analysis of IRAC extended sources, it has been found that light is scattered out 
of the apertures into the array. The effect is most severe in the 5.8 and 8.0 /xm channels. The 
SSC has developed 'best practices' for applying aperture corrections. These corrections are 
dependent upon the IRAC channel and the area within the aperture (or "effective circular 
aperture radii"). All of the reported fluxes (Table 3) have this apperture correction applied, 
which introduces an uncertainty of ~ 10% to the flux densities. No color correction is applied 
to the IRAC data. 

The errors of the IRAC photometry are dominated by two uncertainties: for bright 
sources the uncertainty of the aperture correction is dominant, and for the faint sources the 
uncertainty of the background subtraction becomes important. The errors reported in Table 
3 were determined as the quadratic sum of the two. 

3.2. MIPS 

3.2.1. Observations 

The 24/im observations used the small-field photometry mode. Both the 70/im and 
the 160/im observations used the 3x1 raster map mode, with a 1/8 array step size in the 
cross-scan direction. These provided images of ~ 5' on one side and no less than 2'. 5 on the 
other side, large enough to include some background regions for the sky subtraction, as well 
as to have the central 2'. 5 x 2'. 5 of each map to be fully sampled. For the 3 pairs larger than 
2'. 5, the same observation configuration was applied to each component instead of the total 
pair. For all fields, the same exposure time and same number of cycles were applied, which 
are : (t exp in sees, N cyde ) = (10, 1), (10, 1) and (10, 4) at 24, 70 and 160 /xm, respectively. 
These corresponded to effective integration times of 165.7, 125, and 83.9 seconds for the 
three bands, and yielded average surface brightness rms noise (4<r) of 0.28, 2.0, and 6.0 MJy 
sr _1 in these bands, respectively. 
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3.2.2. Data Reduction 

The BCD images were taken directly from the products of Spitzer Science Center (SSC) 
pipeline (Version 13). These are calibrated images of individual exposures (MIPS Data 
Handbook). The following major artifacts on the BCD images were corrected in our data 
reduction process: 

(1) 'Stripes' in the 70 pim images. These are caused by flat-fielding variations due either 
to the latents of the stim flash or to the long term responsivity variation. All of 
our 70 \im observations were affacted, some more severely than others. Following 
the recommendations in the MIPS Data Manual, we performed self-flat-fielding using 
our own data. This was done by calculating for each detector pixel the 2-cr clipped 
median over all BCD frames within the given observation. Dividing the BCDs by the 
self-flat-field image removed the artificial stripes. 

(2) Dark latents in the 24 fiin images. These are long time scale (up to 10 hours) artifacts 

produced by very bright sources. Two of our observations were affected: AORKEY=14271232 
and AORKEY= 14270720 (see Table 2). In both cases, the bright source which causes 
latents was not in the field of view. The same self-flat-fielding method for the correc- 
tion of the artificial stripes in the 70 /im images was exploited in the correction for 
these artifacts. 

(3) Gradients in the 24 \im images. These non-astronomical gradients are due to poor 
flat fielding. If not corrected, they produce artificial discontinuities in the mosaicked 
images. Most of the 24 /im BCDs were affected by this. We fitted a slope to the sky 
background of each individual BCD and subtracted it. This effectively took out the 
gradients. 

After these corrections, the SSC software MOPEX was used to make the mosaic image out 
of the BCD frames for each 24/im or 70/im MIPS observation. The 160/xm maps were 
taken from the post-BCD products of the same SSC pipeline. Combined mosaic images 
were produced using MOPEX for the 3 large pairs that have MIPS observations of the two 
components separately. 

Flux densities of individual galaxies in pairs were measured from the corresponding 
images. Given the relatively large beams and small separations of pair components, it was 
challenging to carry out these measurements. Two different methods have been exploited: 

(1) Aperture photometry. When only one of the two components in a pair is detectable 
or, if both components are bright in the given MIPS band, the two components are sep- 
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arable, the flux densities were measured using aperture photometry. Polygon apertures 
were used in order to avoid overlapped apertures, with the aperture size changing with 
the galaxy size. The aperture corrections were estimated using two different methods. 
In all cases except for the six 160/im measurements in which the galaxy in question 
is extended in the 160/im band, the aperture corrections were estimated under the 
point source assumption, exploiting the Point Response Function (PRF) presented in 
Appendix B. For the 24/im and 70/xm sources that are significantly larger than a point 
source, this is an underestimation of the real aperture correction, though the fractional 
aperture correction itself decreases with source size. For the extended 160/im sources, 
the apperture corrections were estimated differently. Because of the relatively small 
size of the 160/im images compared to the beam, much of the emission of any galaxy 
that is extended in the 160/im falls out of the 160/im image, and would not be recov- 
ered if an aperture correction based on the point source assumption were applied. In 
order to fully recover the missing flux, the following method was carried out. First, the 
24/xm image (2 times larger than the 160/im image) was convolved with the PRF of 
the 160/im band. Then, under the assumption that the /160//24 ratio is constant across 
the galaxy in question, the desired aperture correction in the 160/im band was esti- 
mated from the ratio between the total flux in the 24/im band and that in the 160/im 
aperture measured on the smoothed 24/xm image. From experiments with individual 
sources, the following aperture correction errors were assigned: 14% for weak sources 
in the 24/xm band, 5% for other 24/xm sources and all 70/im sources, and 10% for all 
160/im sources. 

(2) PRF fitting. This was applied to pairs in which both component galaxies are detectable 
point sources and the two PRFs overlap with each other. In these cases, the flux 
densities of the two components were measured by simultaneously fitting two PRFs 
within a pair, positions and brightness of the two PRFs being the free parameters 
(altogether six of them: xi, yi, peak 1; X2, 1/2, peak 2 ). The average flux measurement 
errors are 14%, 5% and 21% in the 24, 70, and 160/im bands, respectively. 

MIPS data are calibrated assuming a nominal spectrum of a blackbody at 10,000K. This 
is very different from the dust emission in galaxies. Therefore color corrections, calculated 
using the modified blackbody spectra (ft = 2) of temperatures of 100K, 20K and 20K, were 
applied to the 24/im , 70/im and 160/im band flux densities, respectively. The corresponding 
correction factors (multiplicative) are 0.967, 0.923 and 0.954. 

The errors of the MIPS photometry are the quadratic sum of (in the order of impor- 
tance): (1) the aperture-correction/model- fitting error (dominant for all sources except for 
the faintest ones), (2) the uncertainty of the background subtraction, and (3) the rms error 
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of the background. We ignored the photon noise which is never significant. The 4-<r upper- 
limits for the undetected sources were estimated using the quadratic sum of the uncertainty 
of the mean background value and the rms of the background. The apertures addopted in 
these estimates were derived from the K-band Kron radii taken from the 2MASS database. 
However, for a given MIPS band, if the Kron radius of a galaxy is less than a minimum value 
corresponds to the beam, then the aperture is derived using the minimum (8", 20", and 40" 
for 24, 70 and 160/im, respectively). 



4. IR Emission of Paired Galaxies — Images and Catalogs 

Spitzer images at 3.6, 8.0, 24 and 70 fim are presented in Fig.l for all pairs in our sample 
except for KPAIR J1315+6207 (UGC08335a/b). These are gray scale images overlaid by 
contours. The contour levels are Sb g + S x 2 n (n=0, 1, 2, ...), where S& 9 is the local 
background level and S = 0.032, 0.16, 0.28 and 2.0 MJy sr" 1 for the 3.6, 8.0, 24, and 
70/xm band, respectively. The So values are the sample medians of the A-a threshold of the 
corresponding bands. The sizes of all images are 2' x 2' except for the three large pairs: 
J0949+0037, J1453+0317 and J2047+0019, for which a scale bar of length of 2' is given. 
We chose not to present the images in the other 3 Spitzer bands because: (1) the images of 
the 4.5/xm emission have nearly identical morphology as their 3.6/im counterparts; (2) the 
5.8/xm array is the most noisy among the IRAC detector arrays; (3) the angular resolution 
of the 160/im band, FWHM=40", is so coarse that little information can be gained from the 
images in addtion to what is already given by the fi 60 of component galaxies listed in Table 
3. 

The flux densities of KPAIR galaxies, in all 7 Spitzer bands, are listed in Table 3. For 
the 4 IRAC bands, aperture photometry corresponding to physical diameters of 4 and 10 
kpc are also provided. The methods and error estimates are discussed previously in Section 
3. For non-detections, 4<r upper-limits are given. All galaxies in our sample are detected in 
IRAC bands. The detection rates are 47/52, 40/52, and 36/52 for the MIPS 24/im, 70/im 
and 160/im bands, respectively. 

Listed are also two estimates of the total IR (3 — 1000 /im) luminosities. The L TIR is 
derived using the formula of Dale et al. (2005): 

L TIR = 1.559 x L 24 + 0.7686 x L 70 + 1.374 x L 160 (1) 

where L 24 , L 70 , and L 160 are the monochromatic luminosities (vL v ). Another IR luminosity 
is defined using the 24 and 8/im flux densities: 

\og(L IR ) = log(L 24 ) + 0.87(±0.03) + 0.56(±0.09) x log(L 8 /L 24 ) (2) 
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Fig. 1. — Spitzer images of galaxy pairs. The contour levels are Sb g + S x 2 n (n=0, 1, 2, 
...), where Sb g is the local background level and So, the sample median of the 4-cr threshold 
of the corresponding bands, are 0.032, 0.16, 0.28 and 2.0 MJy sr _1 for the 3.6//m, 8.0/im, 
24/im, and 70/im band, respectively. The sizes of all images are 2' x 2' except for the three 
large pairs: J0949+0037, J1453+0317 and J2047+0019, for which a scale bar of length of 2' 
is given. 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 2. — Left: \og{L 2 i/ L TIR ) vs. log(L 8 /L 2 4) plot of spiral galaxies (diamonds: without 
AGNs; eight-point stars: with AGNs) in the KPAIR sample, detected in all Spitzer bands. 
The solid line is a linear regression of the data. The dotted line is the relation found by 
Calzetti et al. (2005) for star formation regions in M51. The dashed line, also taken from 
Calzetti et al. (2005), is the prediction of the model of Dale & Helou (2002) for galaxies. 
Right: \og(L TIR ) vs. \og(L IR ) plot of same galaxies, where \og(L IR ) = log(L 24 ) + 0.87 + 
0.56 x log(L 8 /L 24 ). 
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where L 2 4 and L 8 are again the monochromatic luminosities {vL v ). L IR is an unbiased 
approximation of Ltir- It was derived from the linear regression of the log(L 24 / 'Ltir) 
vs. log(L 8 /L 2 4) correlation using data of 34 spiral galaxies in the KPAIR sample that are 
detected in all Spitzer bands (left panel of Fig. 2), exploiting a similar method that was 
originally developed by Calzetti et al. (2005). The result is nearly identical for a linear 
regression with the galaxies containing AGNs (3 of them) excluded. A comparison between 
log(L TIR ) and log(L IR ) is plotted in the right panel of Fig.2. Indeed there is a very tight 
linear correlation between them, with a dispersion of only 0.11. Given the better detection 
rate of Lir than that of Ltir and the tight linear correlation between the two, we will use 
Lir hereafter whenever total IR luminosities are invoked in calculations. 
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The images and photometry show very diversified IR emission properties among KPAIR 
galaxies. We verify that, except for those with AGN's, paired E galaxies have very low dust 
emission (log(LiR,/L ) ~ 9-10), indicating little star formation occurring in them. Among 
the paired spiral galaxies, the majority have rather moderate IR luminosity (~ 1O 1O L ). 
There are six known AGN's (Table 1), three in S galaxies and three in E galaxies. There 
are four LIRGs (one of them being an AGN), but no ULIRGs, in the sample. Some detailed 
notes on individual pairs can be found in Appendix A. 

5. Star Formation Enhancement 

5.1. Control Sample of Single Late-type Galaxies 

In this section, star formation rates of 39 non-AGN spiral galaxies in the KPAIR sample 
are compared with those of single galaxies in a control sample. In order to have a clean 
comparison, we chose to select a control sample that contains also 39 spiral galaxies, each 
matching a spiral galaxy in the pair sample of the same mass. 

The selection of the control sample was confined to two Spitzer data archives: (1) The 
SWIRE survey of Lockman field and ELAIS-N1 field (covering ~ 20 deg 2 , Lonsdale et al. 

2004) , which provides an IRAC 3.6/im band selected sample of field galaxies (Surace et al. 

2005) . The restriction to the Lockman and ELAIS-N1 fields is because these regions in the 
SWIRE survey have good SDSS spectroscopic coverages (Abazajian et al. 2005); . (2) The 
SINGS survey of nearby galaxies (Kennicutt et al. 2004), which observed a heterogeneously 
selected sample of 75 well-known galaxies, including normal late and early types, AGNs, and 
starbursts. 

The criteria for a galaxy to be considered in the selection of the control sample are: 

(1) has published Spitzer IRAC and MIPS data; 

(2) has spectroscopic redshift, and z < 0.1; 

(3) K s , taken from the K 20 of 2MASS, < 13.5 mag ; 

(4) is a late-type galaxy outside any interacting system; 

(5) does not have a known AGN. 

Paired spiral galaxies with log(M star /M ) > 10.3 were matched with single galaxies in 
the two SWIRE fields. Their Spitzer flux densities were taken from the SWIRE Data Release 



-33- 



2 (Surace et al. 2005). For both IRAC 8fim and MIPS 24/im bands, the Kron fluxes (and 
the associated errors) are adopted. For f 2 4, an additional 15% aperture correction is applied 
(Shupe et al. 2008). SDSS images were inspected and galaxies showing signs of interaction 
or being in pairs were excluded. These single galaxies were separated into two morphological 
types, "S" or "E", using the same scheme as for the classification of paired galaxies (see 
Section 2). Then, each non-AGN spiral galaxy in the KPAIR sample is matched by one of 
the 88 single spiral galaxies so selected, according to the following requirements: (i) The 
single galaxy should have log(M siar ) within ±0.1 of that of the paired galaxy, (ii) Among all 
single spiral galaxies fulfilling the requirement (i), the chosen one should have the minimum 
redshift difference from that of the paired galaxy. It should be pointed out that, despite 
the requirement (ii), there is still a significant difference between redshift distributions of 
the paired galaxies and of single galaxies in the control sample: the medians are z = 0.031 
and z = 0.046, respectively. This is because of the pair selection criterion (1) (see Section 
2) that requires the primaries are brighter K s = 12.5 mag, while the control sample is 
selected from galaxies brighter than K s = 13.5 mag. We argue that the redshift difference 
will not introduce any bias (in particular the Malmquist bias) into the comparisons between 
the two samples that are mass-matched, because all mass-normalized properties (e.g. the 
light-to-mass ratio, SFR/M, etc) shall not depend on redshift in these local samples^, for 
which the cosmic evolutionary effects are negligible. Indeed, the requirement (ii) is arguably 
disposable, and any galaxy fulfilling the requirement (i) could have been included in the 
control sample. In order to assess the uncertainties due to the particular choice of the selected 
control sample, a Monte Carlo analysis (100 repeats) was carried out. In each of the 100 
realizations, an alternative control sample was selected by relaxing the requirement (ii) and 
choosing arbitrarily the match of any KPAIR galaxy among galaxies fulfilling the requirement 
(i). Then the same statistics of the star formation properties that were calculated using 
the official control sample (see the following sections) were repeated using the alternative 
control sample. We confirm that for any of these statistics, the mean of the 100 realizations 
is consistent with the result derived from the official control sample within 1-a. 



2 This insensitivity of mass-normalized properties to any selection effect is the major reason for us to 
choose the one-to-one mass-matching method of control sample selection. 
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There are no single spiral galaxies among the selected SWIRE sources that have log (M/M ) < 
10.3, while 5 KPAIR S galaxies have mass less than this limit. These paired spirals were 
matched by non-interacting normal spiral galaxies in the SINGS sample. There are only 6 
such SINGS galaxies in this mass range. The same requirements that were applied to the 
SWIRE galaxies were applied here, except for NGC 3418 whose mass differs from that of its 
match (J14530523+0319541) by 0.14 dex (the closest match), slightly exceeding the limit of 
0.10 dex. The 8 and 24/im data of SINGS galaxies were taken from Smith et al. (2007). All 
galaxies in the control sample, together with their matches in the KPAIR sample, are listed 
in Table 4. 



5.2. Star Formation Enhancement in Paired Non-AGN Spirals 

The SFR of a galaxy can be estimated from the IR luminosity using the formula of 
Kennicutt (1998) 

SFR (Moyr" 1 ) = 4.510 x 10" 44 x L IR (ergs s" 1 ). (3) 

Note that this estimate is contaminated by the dust emission powered by the radiation of 
old stars. Also, it does not include the UV radiation of young stars that is not absorbed 
by dust. For an average normal spiral galaxy, both biases are at the ~ 30% level (Buat 
& Xu 1996). Under the assumption that they affect the KPAIR sample and the control 
sample in the same way, this estimator is exploited in the comparison between the SFR of 
the two samples. However, when comparing our results with other works using different 
star formation indicators (e.g. the UV continuum and the optical emission lines), the biases 
should be taken into account. 

There are two Lm upper limits in the sample of KPAIR spirals, and two in the control 
sample. In what follows, in any of the statistical analyses involving L/_r, all upper-limits 
were replaced by numbers of half of the values and treated as detections. We confirm that 
there is no significant difference whether these upper limits are or are not included in the 
analyses. 

In Fig.3, a histogram of log(Li R ) (log(SFR)) of non-AGN spirals in the KPAIR sample 
(KPAIR-S) is compared to that of the single spirals in the control sample (CONTROL). 
There is a striking difference between the two distributions: while the distribution of CON- 
TROL has a single prominent peak at log(Li R ) = 9.75 (SFR~ lM yr _1 ), the distribution of 
KPAIR-S has a second peak at much enhanced Li R level of log(Li R ) = 10.65, corresponding 
to SFR~ 8M yr _1 . There is a significant excess of KPAIR S galaxies in the high L IR end. 
Indeed, while 3 non-AGN spiral galaxies in the KPAIR sample are in the LIRG category 
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(log(L IR ) > 11), none of the galaxies in the control sample is a LIRG. The Kolmogorov- 
Smirnov test (K-S test) of the SFR distributions yields a rather low probability of 3.9%, 
or equivalently a rejection at 96.1% confidence level, for the null hypotheses that the two 
samples are drawn from the same population The mean log(LrR) for the KPAIR-S sample is 
log(L m ) = 10.13±0.12, corresponding to a mean log(SFR) = 0.36±0.12. For the CONTROL 
sample, the means are log(L IR ) = 9.84 ± 0.08 and log(SFR) = 0.07 ± 0.08. The Student's 
t-test yields a score of 2.32, corresponding to a probability of 2.6% for the null hypotheses 
that the means of the two samples are equal, consistent with the result of the K-S test. 

Fig. 4 is a comparison of histograms of sSFR, i.e. SFR per mass (SFR/M), of the 
same two samples. It shows a similar shift of the distribution of KPAIR S galaxies toward 
the higher SFR/M bins compared to the control sample, although the difference is slightly 
less prominent than that shown in Fig.3. Indeed, eight KPAIR S galaxies have SFR/M> 
-^Q-iOy r -i w ] 1 ii e on iy one g a i aX y i n the control sample has such high SFR/M. The K-S test 
of the SFR/M distribution finds a low probability for the null hypotheses of 3.9%. The mean 
values of log(SFR/M) are -10.50 ± 0.10 and -10.78 ± 0.08 for the non-AGN spirals in the 
KPAIR sample and in the control sample, respectively. The score of the Student's t-test of 
the means is 2.21, and the probability for the null hypotheses is only 3.3%. 



5.3. Mass Dependence of SFR/M Enhancement in Paired Spirals 



In order to study mass dependence of the SFR and its enhancement in galaxy pairs, we 
binned the spirals both in the KPAIR sample and in the control sample into four mass bins, 
and calculated the means of log(SFR/M) for individual bins. The results are listed in Table 
5 and plotted in Fig.5. Galaxies in the control sample show a clear trend of decreasing sSFR 
with increasing mass, as has already been well documented in the literature (Kauffmann 
et al. 2004; Brinchmann et al. 2004; Schiminovich et al. 2007; Zheng et al. 2007). On 
the other hand, the sSFR of spirals in pairs is nearly constant with mass. In the lowest 
mass bin of 9.7 < log(M/M ) < 10.2, there is no enhancement of the sSFR of the paired 
galaxies compared to that of the control sample. And at the high mass end, in the bin of 
11.3 < log(M/M ) < 11.6, the mean sSFR of the paired galaxies is about 3 times higher than 
that of the control sample. In between, there is a weak enhancement in the two intermediate 
mass bins. 
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Fig. 3. — Histograms of log(Ljft) distributions of the non-AGN spirals in the KPAIR sample 
("KPAIR-S") and of the control sample. Corresponding log(S'Fi?) values are marked on the 
top. 



- 39 - 



14 
12 

10 

8 



2 




KPAIR-S 
CONTROL 



I 



12 -11 -10 

log (SFR/M) (M sun yr 1 ) 



Fig. 4. — Histograms of \og(SFR/M)) distributions of the non-AGN spirals in the KPAIR 
sample ("KPAIR-S") and of the control sample. 



fa 

CO 

o 
Pi 

o 







'Jo 



+ 

CO 

co 



co 
+ 

CO 

co 



O 

O 

U 



CO 

I — I 

< 









CM 








CM 


H 


O 


O 


o 


-H 


-H 


-H 


CM 


b- 


o 


CO 


O 


i> 


O 


H 


o 


i — I 


i — 1 


I— 1 



LO - 


2 


H 


LO* 










o 


LO 




o 


1 — 1 


H 


o 

CM 


co 


o 


CD 




o 


-H 


-H 


23±( 


-H 


i — i 




LO 






LO 


o 


CD 


o 


o 


1 — 1 
1 


H 
i 


H 


H 
1 



lo" 




o 

CM 












1 — 1 


co 




CM 


1 — 1 


CM 


o 

H 


CM 


CD 


o 


o 


o 


-H 


-H 


-H 


HH 


CM 




CD 


co 


oq 


CD 


CD 


CD 


CD 


i — I 


i — 1 


i — 1 


H 


i — I 



LO" 




o 

CM 










o 


oo 


oo 
i — i 


1 — 1 


1 — 1 


i — i 


CM 


o 


o 


o 


o 



i- 1 co 

^ LO 

o o 



o 

CD 
CD 



CM 
CD 



00 CM CD 
CD i— i i— I 



V| VI VI VI 



C2> 

o 



Ol C35 

o o o 



V V V V 

^ 00 CM 

d d h 



-41 - 



-9.5 [ 



-10.0 - 



-10.5 



V 

60 

5 -11.0 



-X: 



-11.5 



KPAIR-S 



CONTROL 



-12.0 L 
9.5 



10.0 



10.5 11.0 

log (M s J (M S J 



11.5 



12.0 



Fig. 5. — Comparison of means of log(SFR/M) of paired and single spirals in mass bins. The 
dotted line marks the mean SFR/M of non-AGN spirals in KPAIR, and the dotted-dashed 
line marks the linear regression of those in the control sample. 
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5.4. Comparison of Non-AGN Spirals in S+S and in S+E Pairs 

Sulentic (1989) and Xu & Sulentic (1991) found significant FIR enhancement in S+E 
pairs in their IRAS studies, under the assumption that the ellipticals in these pairs are FIR 
quiet. However, the ISO observations of Domingue et al. (2003) demonstrated that this 
assumption is invalid. In this section, we address again the question whether spiral galaxies 
in S+E pairs have similar levels of SFR enhancement as those in S+S pairs. 

In Fig.6, the non-AGN spirals in the KPAIR sample is decomposed into two subsamples, 
one for galaxies in S+S pairs (28 galaxies) and the other for galaxies in S+E pairs (11 
galaxies), and the log(SFR/M) distribution of each of them is compared to that of the 
control sample in one of the two panels. Because here the samples being compared do not 
have the same numbers of sources, the distributions are normalized (i.e. in fractions). It 
shows that the distribution of the S galaxies in S+E pairs is not significantly different from 
that of the control sample: the K-S test finds a 91% probability that the two samples are 
drawn from the same population. Therefore, all the enhancement found in Fig.4 is due to 
spirals in the S+S pairs. For the comparison between the spirals in S+S sample and the 
control sample, the K-S test yields a probability of only 2% for the null hypotheses. The 
average values of log(SFR/M) are -10.36 + 0.11 and -10.88 + 0.19 for spirals in S+S sample 
and in S+E sample, respectively. 

In Fig. 7, the mass dependence of SFR/M is plotted for non-AGN spirals in S+S pairs 
and in S+E pairs, separately. No significant enhancement is found for spirals in S+E pairs in 
any mass bin. For those in S+S pairs, they have a similar trend as the total sample (Fig.5), 
showing a rather constant sSFR that is slightly above the mean of the total sample of paired 
spirals. 

Fig.8 is a log(SFR/M)) versus log(M)) plot for individual non-AGN spirals in S+S 
and S+E pairs, compared to their counterparts in the control sample. All spirals in S+S 
pairs are detected. The upper limits for spirals in S+E pairs (two) and in the control sample 
(two) are shown by upside-down triangles and downward arrows, respectively. The two solid 
lines delineate the regions occupied by the LIRGs and ULIRGs. All the four LIRGs are 
in S+S pairs, while there are no ULIRGs in any of the samples. The dotted-dashed line 
marks the result of Brinchmann et al. (2004) for SDSS galaxies. It shows a very similar 
trend in the sSFR vs. mass relation as that revealed by the data of single spirals in our 
control sample, although the SFR in Brinchmann et al. (2004) was derived using the optical 
emission lines data while it was derived using the IR luminosity in this work. The dashed 
line marks SFR/M = l/tnubbie, with the Hubble time tnubbie = 13 Gyr. Galaxies above 
this line have enhanced SFRs compared to a constant SFR over the Hubble time. There is 
indeed an excess of galaxies in S+S pairs above this line, while none of the spirals in S+E 
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pairs has enhanced SFR. 

We define a star formation enhancement indicator, e, for each non-AGN spiral in KPAIR: 



where (SFR/ M)kpair-s and (SFR/M) contro i) are the sSFR of the paired galaxy and that 
of its match in the control sample, respectively. 

Fig. 9 is a plot of means of e in the log(M) bins. For the S+S sub-sample, there is a clear 
mass dependence of e, which can be expressed by its linear regression: 



This relation is confined to the mass range covered by our samples: 10.0 £ log(M/M ) 11.5. 
For the S+E subsample, the mean e is consistent with in all mass bins. 

There is a known dependence of sSFR on local environment in the sense that galaxies in 
higher local density environments tend to have lower sSFR (Kauffmann et al. 2004). Could 
the difference in the SFR/M of spirals in S+S pairs and in S+E pairs be due to different 
mean local densities of the two types of pairs? If S+E pairs are preferentially found in the 
denser environment, then the lower SFR/M of the spirals in these pairs compared to those 
in S+S pairs is just another consequence of the SFR-environment relation. We made the 
following test for this hypothesis. Fig. 10 shows the selection function of the parent sample 
of KPAIR sample (Section 2). In the plot, dots are the 59312 DR3/2MASS galaxies in 
the parent sample, and eight-point stars are the non-AGN spirals in the KPAIR sample. 
The highest redshift of paired spirals is at z ~ 0.06, corresponding to a limiting mass of 
Mii m it = 10 10 ' 9 M . Adopting this limiting mass, around each non-AGN spirals in the 
KPAIR sample we counted neighbors within the parent sample using the following criteria: 
(1) M > M limit ; (2) distance < 2 Mpc; (3) 5V Z < 1000 km s _1 . In Fig. 11, means of the 
neighbor counts (N ne i g hbor) around spirals in S+S pairs and in S+E pairs in the individual 
mass bins are compared with each other. According to Fig. 11, there is no evidence for the 
spirals in S+E pairs having systematically higher local density than those in S+S pairs. 
Therefore the difference in the SFR/M of spirals in the two different types of pairs is not 
due to the local-density dependence. Interestingly, in Fig. 11, there is a clear indication for 
the density/mass correlation, consistent with the literature (Kauffmann et al. 2004). 



Previous studies (Ellison et al. 2008; Woods & Geller 2007) have found that the sec- 
ondaries in minor-merger pairs (mass ratio > 3) have higher SFR enhancement than the 



e = log((SFR/M) K PAiR-s) - log((SFR/M) control ) 



(4) 





5.5. Comparison of Primaries and Secondaries in KPAIR 
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Fig. 8.— The sSFR (SFR/M) versus mass plot for the non-AGN spirals in pairs and 
those in the control sample. The dashed line marks SFR/M = l/t Hubb i e , with the Hubble 
time t Hu bbie = 13 Gyr. The dotted-dashed line is the result of Brinchmann et al. (2004) 
for SDSS galaxies, corrected for the IMF and the Hubble constant differences (the Kroupa 
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primaries. We checked whether this is also true for spirals in KPAIR, which includes only 
major- merger pairs. The answer is negative. As shown in Fig. 12, there is no significant 
difference between the mean SFR/M of primaries and that of secondaries in any mass bins 
studied in this work. The K-S test finds a 62% probability for the null hypothesis that 
SFR/M distributions of the primaries and of the secondaries are drawn from the same pop- 
ulation. 



5.6. Enhancement in One or Two Components? 

In an early IR study of interacting pairs, Joseph et al. (1984) found that among all the 
22 pairs for which they obtained K — L colors for both components, only one component 
of each pair exhibited K — L excess. They interpreted the result as evidence for single- 
component star formation enhancement in interacting pairs. Contrarily, in an IRAS study 
of isolated pairs, Xu & Sulentic (1991) argued that they saw indications of IR enhancement 
in both components in the close interacting (CLO) S+S pairs, though their results were 
not conclusive because most of their CLO pairs are unresolved by IRAS. With the much 
improved angular resolution of Spitzer which resolved all pairs in the KPAIR sample, we 
can address directly the question whether the SFR is enhanced in only one component or in 
both components in close interacting S+S pairs. 

Given our result that only more massive galaxies in S+S pairs have SFR/M enhancement 
(Fig.7 and Fig.9), we picked the 10 pairs (out of 15 in total) whose two components are both 
non-AGN and more massive than lO ia7 M . Fig. 13 is a SFR/M plot of individual pairs, each 
pair in a separate column, showing the log(SFR/M) values of both components as well as of 
their counterparts in the control sample. The pairs are sorted according to the log(SFR/M) 
that is the higher one among the two. 

The results in Fig. 13 can be summarized as the following: 

(1) Massive galaxies (M > lO ia7 M ) in close major-merger S+S pairs have very diversified 

star formation activity levels, from very quiescent (similar to "red and dead" galaxies) 
to strong starburst (e.g. in LIRGs). 

(2) In individual pairs, the SFR/M values of the two components show a certain level of 

concordance ("Holmberg effect"): When one component has a strong star formation 
enhancement, the other is usually enhanced as well (with only one exception: Pair j^l 
= J1043+0645). On the other hand, if one component is a "red and dead" galaxy, the 
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Fig. 12. — Average sSFR's of primaries and secondaries in KPAIR pairs. The mass bins are 
the same as in Fig. 5. The dotted line marks the mean SFR/M of the KPAIR- S sample. 
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Fig. 13.— Log(SFR/M) plot of 10 S+S pairs with M > lO iar M . Filled circles represent 
paired galaxies, and open squares represent their counterparts in the control sample. Pairs 
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other one usually shows no sign of star formation enhancement, eitheiQ 

Fig. 14 confirms the "Holmberg effect" between the SFR/M of the two pair components. 
While the values of log(SFR/M) of the two components are highly correlated, no correlation 
is found between log(SFR/M) values of their matches in the control sample, as is expected. 

The result (2) is apparently in contradiction with that of Joseph et al. (1984). This 
might be explained by the differences between our study and that of Joseph et al. (1984). 
Several pairs in the sample of Joseph et al. (1984) belong to the category of minor-mergers 
(e.g. Arp 283 and Arp 294), for which indeed only one component (the secondary) is usually 
enhanced (Woods & Geller 2007). Furthermore, the study of Joseph et al. (1984) was con- 
fined to the nuclei of the interacting galaxies, and the K— L excess is only sensitive to very hot 
dust emission from compact starbursts. It is possible that these compact nuclear starbursts 
have shorter time scales, and therefore significantly lower chance to occur simultaneously in 
both components, compared to the star formation enhancement over entire galaxy bodies as 
probed by our Spitzer observations. 

Fig. 15 is a plot of log(SFR/M) versus neighbor counts (N neighbor , as defined in Section 
5.4) for the same paired galaxies plotted in Fig. 13. There is no discernible dependence of 
log(SFR/M) on N nei g hbor in the plot. Therefore, it is something other than the local density 
that determines whether the two galaxies in an S+S pair should or should not have enhanced 
star formation activity. 



5.7. Separation and SFR/M Enhancement 

Previous studies (Xu & Sulentic 1991; Barton et al. 2000; Lambas et al. 2003; Nikolic 
et al. 2004; Alonso et al. 2004; Woods et al. 2006; Barton et al. 2007; Ellison et al. 2008) 
found that the star formation enhancement in pairs of separation ^ 20 h _1 kpc is much 
stronger than those of larger separations. Using our sample, we can address the question 
whether there is still dependence of the SFR on separation for pairs within this separation 
limit. We define a normalized separation parameter: 

SEP = — ^ (6) 
ri + r 2 

where s is the projected separation, and i\ and r 2 are the K-band Kron radii (taken from 
2MASS) of the primary and the secondary, respectively, in the same units as those of s (kpc 



3 This is in agreement with the lack of SFR/M enhancement for spirals in S+E pairs. 
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Fig. 14. — Correlation plot of Log(SFR/M) of two components in the S+S pairs in Fig. 13. 
Filled circles represent the pairs, and open squares represent the their counterparts in the 
control sample. 
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in Fig. 13. 
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or arcsec). In the ideal case of a pair of two round galaxies, the two component galaxies 
overlap with each other when SEP < 1. 

Fig. 16 is a histogram of SEP distribution of non-AGN spirals in KPAIR. The median 
is at SEP ~ 1 (mean SEP = 1.05 ± 0.40). Fig.17 is a plot of mean log(SFR/M) v.s. log(M) 
of spirals in S+S pairs, separated into two subsamples of SEP < 1 and SEP > I. For both 
subsamples, the log(SFR/M) v.s. log(M) relation scatters around the mean of the total 
sample without any obvious trend, and no significant difference between the two subsamples 
is detected. This seems to suggest that the separation is not an important parameter any 
more once the two galaxies are close enough. There might be several conflicting factors 
affecting the star formation activity vs. SEP relation in close major-merger pairs. On the 
one hand, colliding pairs with SEP < 1 may undergo collisionally triggered starbursts in the 
regions where the two galaxies overlap, as in the case of the Antennae Galaxies (Xu et al. 
2001). On the other hand, Gao & Solomon (1999) found a correlation between molecular gas 
content and the pair separation, suggesting a progressive gas depletion due to prolonged star 
formation activity. It should also be noticed that SEP is derived from the projection of the 
real, 3-dimensional separation. Any dependence on the true separation can be significantly 
disturbed by the projection effect. 

5.8. Contribution of KPAIR galaxies to the Cosmic SFR Density 

In this subsection we estimate the contribution of galaxies in close major-merger pairs, 
as defined by the pair selection criteria in Section 2, to the total cosmic SFR density in the 
local universe. Assuming the contribution for E galaxies is negligible, this can be estimated 
as follows: 



where ip(M) is the mass function of KPAIR galaxies and f s the S fraction (N s /N). Because 
the mass dependence of (SFR/M) K pair-s is rather flat, we assume it is constant and equal 
to the mean SFR/M of the KPAIR-S sample: (SFR/M) KP air-s = lO -10 - 50 yr -1 . The mass 
function ip and the S fraction f s are taken from Domingue et al. (2009). The result of the 
integration is Pkpair = 2-54 x 10~ 4 M yr _1 Mpc -3 , which is 1.7% of total cosmic SFR 
density in the local universe (p' = 0.015 M yr -1 Mpc -3 ; Yun et al. 2001). 

It should be pointed out that, because of the criterion on pair separation (r > 5h -1 
kpc), mergers already coalesced are missing in the KPAIR sample. This population includes 
the majority of ULIRGs and many LIRGs. Assuming that all ULIRGs and 70% LIRGs are 
coalesced mergers, and estimating their densities using the IR luminosity function of Yun et 
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Fig. 16. — Histogram of the SEP distribution of non-AGN spirals in the KPAIR sample. 
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Fig. 17. — Average sSFR's of non-AGN spirals in S+S pairs, separated into two subsamples 
of SEP < 1 and SEP > 1, respectively. The mass bins are the same as in Fig.5. The dotted 
line marks the mean SFR/M of the non-AGN spirals in the KPAIR sample. 
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al. (2001), we found that the total SFR in these sources is p = 2.77 x 10" 4 (M yr" 1 Mpc" 3 ), 
i.e. nearly the same as that in the KPAIR galaxies. Therefore, in the z=0 universe, the total 
SFR in close major mergers is p = 5.31 x 10~ 4 (M yr _1 Mpc~ 3 ). This is only 3.5% of the 
local cosmic SFR density, a truly negligible contribution indeed. 

6. Discussion 

6.1. Dependence of SFR Enhancement on Companion's Morphological Type 

The star formation enhancement found in interacting galaxies is often explained in 
terms of gas inflow caused by gravitational torques of interaction-induced bars (Hernquist 
& Barnes 1991; Barnes & Hernquist 1996). However, this theory cannot explain our result 
that the star formation enhancement depends on the morphological type of the companion 
galaxy: while spirals in S+S pairs show significant enhancement, those in S+E pairs have 
star formation activity comparable to that of single spirals. It appears that, in addition to 
pure gravitational effects, some other factors related to the companion must play important 
roles in the interaction induced star formation. 

A related result was reported by Park & Choi (2009) in a study of the dependence of 
physical parameters of SDSS galaxies on small scale and large scale environments. These 
authors found that the SFR of late-type galaxies is enhanced when the nearest neighbor 
is also a late-type, but reduced when the neighbor is an early-type. They suggested that 
the hot gas halo of an early-type companion can suppress the SFR of a late-type galaxy, 
through hydrodynamic effects such as ram pressure stripping, viscous stripping and thermal 
evaporation, analogous to what is being encountered by late-type galaxies in clusters (Boselli 
& Gavazzi 2006). This interpretation can be applied to our result of low sSFR enhancement 
of spirals in the S+E pairs. Verdes-Montenegro et al. (2001) argued that a similar mechanism 
might be responsible for the depressed SFR of galaxies in compact groups (but see Rasmussen 
et al. 2008). The X-ray observations of four S+E pairs by Gruntzbauch et al. (2007) indeed 
showed evidence for extended X-ray halos in the E components, lending support to this 
hypothesis. 

It should be noticed, however, that our result of no SFR enhancement for spirals in S+E 
pairs is based on a small sample of 11 S+E pairs, and therefore should be confirmed by future 
studies using larger samples. Furthermore, there are known exceptions of strong starbursts 
in S+E pairs in the literature, such as NGC 3561 A (a LIRG) in Arp 105 (Due et al. 1997). 
It will be worthwhile to investigate why galaxies like NGC 3561 A behave differently. One 
noticeable feature of NGC 3561A is the segregation of atomic and molecular gas: the HI 
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gas is completely displaced out of the disk, and the nuclear starburst is supported by pure 
molecular gas (Due et al. 1997). Another spiral galaxy in an S+E pair that has a similar 
atomic-molecular gas segregation is NGC 1144 in Arp 118 ("Yin- Yang Galaxy", Appleton et 
al. 2003), also a LIRG, though in this case the nucleus is an AGN rather than a starburst. 

6.2. "Holmberg Effect" 

The "Holmberg effect" on SFR/M of massive S+S pairs (Fig. 14) is in agreement with 
the result of Kennicutt et al. (1987) derived from the integrated H a fluxes, and that of 
Hernandez- Toledo et al. (2001) based on the (B-V) colors. Apparently, the effect is present 
in star formation indicators of very different time scales (~ 10 7 yrs for the H a emission, 
~ 10 8 yrs for the IR emission, and ~ 10 9 yrs for the (B-V) color). On the other hand, 
it has been found only in global star formation indicators for entire galaxies, but not in 
those for the nuclear star formation activity (Joseph et al. 1984). Interestingly, the SFR 
dependence on interaction parameters has been invoked to explain both the presence of the 
correlation between the SFR of the two components (Kennicutt e al. 1987), and the absence 
of it (Joseph et al. 1984). However, we have shown that for spirals in close major-merger 
pairs, interaction parameters such as the separation are not important factors in determining 
whether a galaxy has enhanced SFR/M or not (Fig. 17). It is possible that the concordant 
star formation behavior of galaxies in a close major- merger pair is dictated by the local 
environment within/around the dark matter halo (DMH) surrounding the pair. But, as 
shown in Fig. 15, the level of star formation activity of these pairs depends very little on the 
local density. It is possible that the SFR is suppressed in those quiescent S+S pairs because 
there is diffuse hot IGM gas in the DMH, in a similar way as what may be happening in the 
S+E pairs (see Section 6.1). Or, in a related scenario, it might be because the DMH's of 
these quiescent S+S pairs have no "cold streams" of IGM gas (Keres et al. 2009) to fuel the 
star formation in the component galaxies. It will be worthwhile to confirm or refute these 
speculations in future studies. 

6.3. Dependence of SFR Enhancement on Mass 

Our result indicating a lack of SFR enhancement in low mass late-type interacting 
galaxies is in agreement with observations of Brosch et al. (2004) and Telles & Maddox 
(2000). In the theory proposed by Mihos et al. (1997), this is due to the fact that these 
galaxies do not have sufficient disk self-gravity to amplify dynamical instabilities, and this 
disk stability in turn inhibits interaction-driven gas inflow and starburst activity. 



- 61 - 



On the other hand, studies including minor-mergers (Woods & Geller 2007; Ellison et 
al. 2008; Li et al. 2008) found significant SFR enhancement in low mass (M < 1O 1O M ) 
interacting galaxies. It appears that, in a minor-merger pair, a low mass galaxy can have 
much stronger SFR enhancement when its companion is much more massive (Woods & Geller 
2007). 

6.4. Overall SFR Enhancement in Major-Merger Pairs 

We have found that, for close major mergers, only massive (M ^ 1O 1O ' 5 M ) galaxies in 
S+S pairs have significant star formation enhancement. These galaxies are less than 30% in a 
K-band selected sample of close major mergers (Domingue et al. 2009), and some of them are 
locked in low sSFR pairs (Fig. 13). Therefore, even for spiral galaxies in close major-merger 
pairs which harbor most merger-induced star bursts in the universe, the star formation 
enhancement due to galaxy-galaxy interaction is still confined to a small sub-population. 
This is consistent with the observations of Bergvall et al. (2003) and the simulations of Di 
Matteo et al. (2008), both argued that mergers are not very efficient in triggering significantly 
enhanced star formation. This is also consistent with the low contribution of major-merger 
galaxies to the cosmic SFR density in the local universe (Section 5.8) and in the uninverse 
of intermediate redshift (z~ 0.24 — 0.80, Jogee et al. 2009). 

7. Summary 

We present Spitzer observations for a sample of close major-merger pairs of galaxies, 
selected from a 2MASS/SDSS-DR3 cross-match. The scientific goals are (1) studying the 
star formation activity in these galaxies and (2) setting a local bench mark for the cosmic 
evolution of close major mergers. The Spitzer KPAIR sample (27 pairs, 54 galaxies) includes 
all spectroscopically confirmed spiral-spiral (S+S) pairs and spiral-elliptical (S+E) pairs 
in a parent sample that is complete for primaries brighter than K=12.5 mag, projected 
separations of 5 < s < 20h~ 1 kpc, and mass ratios < 2.5. There are 42 spiral galaxies and 
12 elliptical galaxies in the sample. These galaxies harbor 6 known AGNs, 3 in spirals and 
3 in ellipticals. 

Spitzer observations include images in the four IRAC bands at 3.6, 4.5, 5.8 and 8.0 
/im, and the three MIPS bands at 24, 70 and 160 fim. They show very diversified IR 
emission properties among KPAIR galaxies. Among the paired spirals, the majority have 
rather moderate IR luminosity (~ 1O 1O L ). There are four LIRGs (~ 10% of KPAIR-S 
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subsample), but no ULIRGs. The SFR, estimated using the IR luminosity, and the sSFR 
of non-AGN spirals (39 of them) in KPAIR are compared to those of single spirals in a 
control sample. Each of the 39 galaxies in the control sample matches a non-AGN spiral in 
KPAIR with the same mass (estimated from the K-band luminosity). The following results 
are found: 

(1) The mean SFR of non-AGN spirals in the KPAIR sample (KPAIR-S) is significantly 
enhanced compared to that of the single spirals in the control sample. The means 
of log(SFR) of the KPAIR-S galaxies and of the galaxies in the control sample are 
0.36±0.12 and 0.07±0.08, respectively. And the means of log(SFR/M) of the KPAIR-S 
galaxies and of the galaxies in the control sample are —10.50 ±0.10 and —10.78 ±0.08, 
respectively. The K-S test rejects at the 96.1% confidence level the null hypotheses 
that the two samples are drawn from the same population. 

(2) When separating the non-AGN paired spirals into those in S±S pairs (28) and in S±E 

pairs (11), only the former show SFR/M enhancement whereas the latter do not. 
The means of log(SFR/M) of the spirals in S±S pairs and of those in S±E pairs are 
-10.36 ± 0.11 and -10.88 ± 0.19, respectively. 

(3) The SFR/M enhancement of spirals in S±S pairs is highly mass-dependent: only those 

with M > lO ia5 M show significant enhancement, whereas relatively low mass (M ~ 
1O 1O M ) spirals in S±S pairs have about the same SFR/M as their counterparts in the 
control sample. 

(4) We define e as the SFR enhancement parameter, e = log((SFR/M) K pAiR-s)— log((SFR/M) contro i), 

where (SFR/M) K pair-s and (SFR/M) contro i are the SFR/M of a paired spiral and 
that of its match in the control sample, respectively. For spirals in the S±S sub- 
sample, there is a strong linear dependence of e on log(M), specified as < e >s+s= 
0.03(±0.14) + 0.47(±0.15) x (log(M/10 10 M o )). The relation is valid for the mass range 
of 10.0 < log(M/M ) < 11.5. 

(5) For spirals in the KPAIR sample, which includes only close major mergers, there is no 

systematic difference between the log(SFR/M) of spirals with SEP < 1 and those with 
SEP > 1, SEP being the normalized separation parameter (SEP = s/(ri +r2)). Also, 
there is no significant difference between the means of log(SFR/M) for the primaries 
and the secondaries. 

(6) There is evidence for a correlation between the global star formation activities (but not 

the nuclear activities) of the component galaxies in massive S±S major-merger pairs 
(the "Holmberg effect"). 
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(7) The contribution of KPAIR galaxies to the cosmic SFR in the local universe is P^pair = 
2.54 x 10~ 3 M yr" 1 Mpc~ 3 . This is 1.7% of total cosmic SFR density in the local 
universe. Adding the SFR in mergers already coalesced, which are missed by the 
KPAIR sample and may include many ULIRGs and LIRGs, the total SFR in close 
major mergers is p' = 5.31 x 10~ 3 M yr _1 Mpc~ 3 . This is only 3.5% of the local 
cosmic SFR density. 
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APPENDIX 



A. Notes on Individual Pairs 

J0020+0049 The FIR emission of this S+E pair is dominated by the S component. The 
E component was only barely detected in the MIPS 24 fxm band, and undetected in the 70 
and 160 \xm bands. 

J0109+0020 Neither galaxy in this S+E pair is detected in any of the MIPS bands, a unique 
case in the sample. The S component shows an arm-like feature in the optical images, but 
was classified as an E galaxy by the automatic classification routine based on the optical 
color and light concentration. Apparently it is a 'red and dead' galaxy, perhaps of early S 
type. 

J0118-0013 Both components of the S+S pair are well detected by IRAC and MIPS. The 
western component has a narrow line AGN (Hao et al. 2005). In its optical image there 
is a blue, jet-like feature pointing to the companion galaxy. It is a "luminous IR galaxy" 
(LIRG), with log Ltir/Lq = 11.41, dominating the total dust emission of the pair. From 
aperture photometry of the IRAC bands, the AGN contributes ~ 40% of the dust emission 
of the galaxy. The eastern component looks like a normal late-type spiral, contributing only 
~ 12% of the Ltir of the pair. 

J021 1-0039 This S+S pair is KPG 058 (Karachentsev 1972). Both galaxies are seen edge- 
on. The western component (an Sbc galaxy) has a narrow line AGN (Hao et al. 2005). It 
dominates the dust emission of the pair. 

J0906+5144 Another KPG pair, KPG 185. The western component is classified as E type 
in our scheme, though according to NED it is an Sa galaxy. The eastern component is an 
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S type with a narrow-line AGN (Hao et al. 2005). The IRAC 8 /im band image shows a 
nucleus+ring morphology, with most of the emission from the ring. 

J0937+0245 This S+E pair is Arp 142 (= VV 316), classified as 'ring galaxy' by RC2 (de 
Vaucouleurs et al. 1976). The S component (NGC 2936) is very disrupted. The optical/NIR 
morphology looks like a bird head, with the nucleus being the eye. The tidal tails, which form 
the 'neck' of the bird, contain several star forming regions bright in the IR (see the 8 and 24 
fim image). However, most of the dust emission is still confined within the nucleus+disk (i.e. 
the bird head) region. The E component (NGC 2937) is undetected in the 70 and 160/im 
bands. 

J0949+0037 This S+S pair is KPG 216. It is one of the three pairs in our sample larger 
than 2'. 5. The two relatively low mass galaxies (~ 1O 1O M ) are well separated from each 
other. The Spitzer observations were carried out for the two components separately, and the 
final maps are coadds of these separate observations. Both components are well detected by 
IRAC and MIPS. And both show extended IR emission through out the entire discs. There 
is a strong outer disc starburst on the east side of the eastern component (NGC 3023), whose 
FIR luminosity is comparable to that of the nucleus. 

J1020+4831 The eastern component (E type) is a strong radio source (4C +48.29), clas- 
sified as an AGN in the literature according to NED. It is a rather weak FIR source, only 
marginally detected in the 24 \im band, and undetected in the 70 and 160 fim bands. The 
western component (S type), detected in all 7 Spitzer bands, dominates the dust emission of 
the pair. 

J1027+0114 This close S+E pair is actually in a triplet (Karachentseva et al. 1988). The 
third galaxy, west of the pair, is about 2' away from the pair center. The southern component 
(S type) is a strong IR source. The northern component (E type) is detected marginally by 
MIPS in the 24 fim band. 

J1043+0645 Both components of this S+S pair are detected by IRAC and MIPS. The two 
components have nearly equal mass (the mass ratio derived from the K band lumiosities is 
Mi/M 2 = 1.2), but rather uneven dust emissions. The L T ir of the western component is 
more than 5 times of that of the eastern component. 

J1051+5101 This is KPG 253, which appears to be a very close S+E pair. The projected 
separation of the two components (s = 4.7 h _1 kpc) is actually less than the 5 hr l kpc 
lower boundary of the separation criterion. However, since the separation is only 6% off 
the boundary, we chose to keep the pair in this work. It is in the center of a cluster. Both 
galaxies are massive (> lO n M ). The westhern galaxy is classified as E while the eastern 
galaxy as S. At low surface brightness levels, a ring like structure shows up around the 
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nucleus of the S component both in the optical and in the 8/xm images. The pair is detected 
by IRAS with f 60 = 0.78 Jy. Interestingly, both the 24/wn and 70/zm band emission contours 
coincide with the nucleus of the E component while there is little emission detected in the 
nuclear region of the S component. This indicates that most of the star formation in the 
pair is occuring in the nucleus of the more massive E component, perhaps due to an ISM 
transfer from the S component. The total SFR derived from the IR emission is quite low, at 
the level of ~ lM G /yr, very different from IR selected mergers. No AGN or any detectable 
radio source has been found in the pair (van Driel et al. 2000). 

J1202+5342 The weastern component (S type) dominates the dust emission the S+E pair. 

J1308+0422 =UGC 8217. Both galaxies in this S+S pair are well detected in all 7 Spitzer 
bands with nearly equal L T ir 4 1O 9 L ). 

J1332-0301 Both spiral galaxies are detected by IRAC and MIPS. The western component, 
with the lower mass and more compact morphololgy among the two, has warmer (i.e. higher 
f 70 /fi6o ratio) and stronger dust emission than the eastern component. 

J1346-0325 The two galaxies of this S+E pair are well separated from each other. Both 
are detected by IRAC and MIPS, with comparable, moderate dust emission. The eastern 
component (an E type) has a Sy2 nucleus (Maia et al. 2003). 

J1400+4251 This S+S pair is a bright IRAS source (f 60 = 2.32 Jy). The combined L Tm 
(= 1.3 x 10 n L o ) makes the pair a LIRG. The two spiral galaxis, both detected by IRAC 
and MIPS, are about equally bright in the dust emission. When separated, each is slightly 
fainter than a LIRG. 

J1425+0313 This S+E pair is VIII Zw 415. The western component (E type), a moderate 
IR source, contains a broad line AGN (Hao et al. 2005). 

J1433+4004 =KPG 426. These are two IR bright spiral galaxies, both are detected by 
IRAC and MIPS. There are bright hot-spots in the 8 /im image of the southern galaxy. As 
a close-interacting S+S pair, it has been studied by Xu et al. (2001) using ISOCAM. 

J1453+0317 One of the three pairs in our sample that are larger than 2'. 5. These are two 
well separated S-type galaxies, both are detected by IRAC and MIPS. The eastern galaxy 
displays a ring formation in both 8 and 24/im images. The western galaxy was detected in 
the ISOPHOT 170 micron Serendipity Survey (Stickel et al. 2004), with f 170 = 1.94 ± 0.58 
Jy. This is significantly lower than the MIPS 160 fxm flux: fi 60 = 3.03 ± 0.32 Jy. 

J1506+0346 An S+S pair. The western component (IC 1087), classified as SO-a (via NED), 
is a weak IR source. Most of the dust emission is from the eastern component (UGC 09710). 
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J1510+5810 These are two very close S-types. The source about 1' west of the pair is a 
star. Most of the dust emission is due to the eastern component. It was detected in the 
ISOPHOT 170 micron Serendipity Survey (Stickel et al. 2004), with f 170 = 0.87 ± 0.26 Jy. 
This is slightly lower than the MIPS 160 //m flux: fi 60 = 1.16 ± 0.12 Jy. 

J1528+4255 =1 Zw 113. These two spiral galaxies show signs of interaction. Most of the 
dust emission is found in the western component (NGC 5934). 

J1556+4757 The western component (S type) dominates the dust emission in the S+E 
pair. The source between the two galaxies is a star. 

1602+4111 =KPG 479. Both spiral galaxies are IR bright. They show signs of interactions. 

1704+3448 These are two closely-interacting, star-forming spiral galaxies. The northern 
component is a LIRG. 

J2047+0019 =KPG 548. The S+E pair is one of the 3 pairs in the KPAIR sample with 
size > 2'. 5. The western component (S type), well resolved in both the IRAC and the MIPS 
bands, shows a nucleus+ring structure. In the MIPS bands, the ring is rather faint. The 
eastern component (E type) is only marginally detected in the 24/im band, undetected in the 
70 and 160 /im band. This pair has been studied by Domingue et al. (2003) with ISOPHOT. 

J1315+6207 =UGC08335=Arp 238=KPG 369=VV 250. This S+S pair is not displayed in 
Fig.l. It was observed successfully by Mazzarella et al. (2009) in the IRAC bands, but the 
MIPS observations failed. The IRAC images (Mazzarella et al. 2009) and the KAO image at 
100/im (Bushouse et al. 1998) show that the eastern component dominates the IR emission 
of this LIRG pair. The west component was undetected in the KAO observation (Bushouse 
et al. 1998). The radio continuum observations at 4.85 GHz (Condon et al. 1991) revealed 
a flux ratio between the two components of 4.4 (22mJy/5mJy), close to the 8/xm flux ratio 
of 3.5 (187mJy/54mJy) measured by Mazzarella et al. (2009). Compared to the IRAS data, 
Condon et al. (1991) found a FIR/radio ratio index q = 2.74 for the pair, very close to the 
mean (< q >= 2.64 ± 0.16) for star forming galaxies. Assuming that both components have 
the same FIR/radio ratio (e.g. Hattori et al. 2004), we divided the total IR luminosity of 
the pair, \og(L T iR/ L Q ) = 11.74 (Mazzarella et al. 2009), and found \og(L TIR / L Q ) = 11.65 
for the east component and \og(L T m/ L Q ) = 11.01 for the west component. Therefore, both 
galaxies qualify as LIRGs. 
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B. MIPS Point Response Functions (PRFs) 



The 24/j.m PRF was taken from Engelbracht et al. (2007). It was generated by smooth- 
ing the standard Spitzer TinyTim PSF (Krist 1993) with a 4". 41 ('1.8-pixel') boxcar. Simi- 
larly, the 70 [im PRF was taken from Gordon et al. (2007), which is a TinyTim PSF smoothed 
by a 13". 30 ('1. 35-pixel') boxcar. As shown in Fig.A-1, good agreements were found between 
our data and the above model PRFs. 

On the other hand, the comparison between our data in the 160/im band and the 
standard PRF ("STinyTim+ 1.6-pixel boxcar") of the MIPS team (Standsberry et al. 2007) 
showed an excess in the first Airy-ring in our data (Fig.A-2). In order to better fit the data, 
we adopted an empirical PRF. It is the standard PRF plus a ring, which has a truncated 
Gaussian profile: 



S 




i 



(r < r 2 ) 
(r > r 2 ); 



(Bl) 
(B2) 



= 



where Si = 0.05, r 1 = 5 pixel, r 2 = 8 pixel and a x = 4 pixel. The pixel size is 8" for the 
160/im map. The comparisons between 160/im PRFs and the data can be found in Fig.A-2. 




Fig. A-l. — Left The 24/im point response functions (PRFs). right: The 70 jim PRFs. 
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Fig. A-2.— The 160/xm PRFs. 



